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a^ , 

Considei|3 a weig ht p{dz) := e-^(^)dz on an interval F = [A,B] C R, with rational 
C . logarithmic derivative and subjected to the following boundary conditions: 

in . V' = ^= ^1 ' lim f(z)e-^^'h'' = for all k > 0. (1) 

On the ensemble 

Q>^ . 'Hn = {N X Hermitean matrices} 

o ■ 

Q ' define the probability 

^1 P{M edM)=CNe-^'^^^'^dM, (2) 

^ ! with the Haar measure dM on TYat; the latter can be decomposed into a spectral part 

(radial part) and an angular part: 

N 

dM ■=\{dMii W (dReMijdlmMij) = A'^{z)dzi...dzN dU, (3) 
where A(2;) = H (-^i ~ -^j) is the Vandermonde determinant. This lecture deals 



^> , l<i<j<7V 

rS \ with the following probability for E C F: 



P(spectrum M G E, with M e H^) = {^"^^iiSlf^ 

JpnA^{z)IIi p{dzk) 

As the reader can find out from the excellent book by Mehta [], it is well known 
that, if the probability P{M G dM) satisfies the following two requirements: (i) 
invariance under conjugation by unitary transformations M i-^ UMU~^, (ii) the ran- 
dom variables Ma, Re Mij, Im Mij, I < i < j < N are independent, then V{z) is 
quadratic (Gaussian ensemble). For this ensemble and for very large A^, the probabil- 
ity P(an eigenvalue G dz) tends to Wigner's semi-circle distribution on the interval 

^Appeared in: Group 21, Physical applications and Mathematical aspects of Geometry, Groups 
and Algebras, Vol.11, 835-852, Eds.:H.-D. Doebner, W. Scherer, C. Schulte, World scientific, Singa- 
pore, 1997. 
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[—\/2N, \/2N]. Upon normalizing the distribution such that the average spacing 
between eigenvalues equals one, one finds different limits for ^00, according to 
whether one considers the part of the spectrum near z — Q (bulk scaling limit) or 
near the edge \/2N (edge scaling limit). The first situation leads to the Sine kernel 
and the second one to the Airy kernel. 

It is also known that the following probability can be expressed as a Fredholm 
determinant, for both, N finite and N infinite. A very sketchy proof of this fact will 
be exposed in this lecture. Therefore the real issue is to compute for £■ C F, 



P (exactly k eigenvalues in E) 
for the kernel 



(-21 



(5) 



K''{y,z) = K{y,z)lE{z). (6) 

Below are a few important examples of kernels, the first one being one side of the 
Christoffel-Darboux formula for classical orthogonal polynomials: 



e 2 



i(y(2/)+y(z)) 



N 



^Pk-i{y)Pk~i{z) 

k=l 



2 Jo 



y- z 



y + z 



rx 

/ A(x + y)A(x + z)dx 
Jo 



("Christoffel- 
Darboux" kernel) 
(Sine kernel) 

(Airy kernel) 



xJ„{x^Jy)Ji,{x\fz)dx 



(Bessel kernel). 

(7) 

As a feature, common to all of them, note that, for instance, the orthogonal polyno- 
mials Pk{y), the exponential function e~*^^ ie*"^^, the Airy function A{x + y) and the 
Bessel function Jy{x^) are all eigenfunctions of second order problems. 

Random matrices provide a model for excitation spectra of heavy nuclei at high 
excitations (Wigner [], Dyson [] and Mehta []). Since the energy levels are gov- 
erned by the spectrum of a quantum mechanical Hamiltonian on a Hilbert space, 
it is reasonable upon truncation to assume that the Hamiltonian be represented by 
large matrices and, from a statistical mechanical point of view, by an ensemble of 
matrices, possibly satisfying some symmetry properties. In their analysis of nuclear 
experimental data. Porter and Rosenzweig [] observed that the occurrence of two 
levels, close to each other, is a rare event (level repulsion) , showing that the spacing 
is not Poissonian, as one might expect from a naive point of view; this lead Wigner 
to his so-called surmise. 

In their pioneering work, Jimbo, Miwa, Mori and Sato [] have shown some (p, q)- 
variables derived from the sine kernel satisfy a certain Neumann-like completely inte- 
grable finite-dimensional Hamiltonian system. Tracy and Widom [] have successfully 
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used functional-theoretical tools to compute the level spacing distributions for a more 
general class of kernels, always yielding Neumann-type systems. In the case where 
is a semi-interval, they find a distribution whose density satisfies a Painleve type 
equation. Deift, Its and Zhou [] have used the Riemann-Hilbert approach to find 
the precise asymptotics for the distributions above. Random matrices have come up 
in the context of statistical mechanics and quantum gravity; see []. Random matrix 
ideas play an increasingly prominent role in mathematics: not only have they come 
up in the spacings of the zeroes of the Riemann zeta function, but their relevance has 
been observed in the chaotic Sinai billiard and, more generally, in chaotic geodesic 
flows; Sarnak [] conjectures that chaos leads to the "spectral rigidity", typical of 
the spectral distributions of random matrices, whereas the spectrum of an integrable 
system is random (Poisson)! 

The main question of this lecture is the following: What is the connection of 
random matrices with integrable systems ? Is this connection really useful ? Re- 
member an integrable system is a time evolution involving commuting vector fields, 
parametrized by times ti,t2, ... . Celebrated examples are the KP hierarchy and the 
Toda lattice. 

The approach in this lecture, based on joint work with M. Adler and T. Shiota 

[], is novel and different from the previous studies, mentioned above. Here we get a 
system of partial differential equations directly for the probability (5) with k = 0, 
rather than a Neumann system for some complicated auxiliary variables. Moreover 
our PDE's are nothing else but the KP hierarchy for which the i-partials, viewed as 
commuting operators, are replaced by non-commuting operators in the endpoints Ai 
of the interval E under consideration. When the boundary of E consists of one point, 
one recovers the Painleve equations, found in [], at least for some of the kernels (7); 
from our work, it also appears that some of the Painleve equations can be viewed as 
the KP equation in non-commutative operators. 

Invoking a different integrable system, H. Peng [] has extended these methods 
to symmetric ensembles, instead of Hermitean ensembles. Symplcctic ensembles will 
require the use of still a different integrable system. Picking the 2-Toda lattice leads 
to equations for the distribution of the spectrum associated with two coupled random 
matrices; see []. 

The first half of this lecture discusses the finite Hermitean matrix ensembles, which 
relate to discrete integrable models, while the second deals with infinite ensembles, 
which arise in the context of continuous systems. The first part serves, to a large 
extent, to understand and motivate the methods developped in the second part. 

1 Finite ensembles at the crossroads of KP, Toda and Vira- 
soro 
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1.1 Introducing time in the probability 

On the ensemble we consider probability (2), in which we introduce "time", 
seemingly as an artifact: 

P(M„ e dM) = c„(i)e-^^W+Sr*'^^^dM 

and thus 

Jpn A2(z) nfe=i Pt{dZk) Tn[t) 

where 

1.2 Tn{E,t) satisfies the KP hierarchy 

Introduce orthonormal and (monic) orthogonal polynomials, Pi{z) := pi{E,t,z) and 
Pi{z) := Pi{E, t, z) respectively, with regard to the density pt, on E G F C R, 



{Pi,Pj)E,t^ Sij, {pi,pj)E,t^ hiSij, Pi^yhiPi. (10) 
The matrix m„ of moments fiij{E,t), at time t 

mn{E, t) := {p,,{E, i))o<,,<„_i := {{z\ ^^■)^,*) , (H) 

can be expanded in Schur polynomials Sk{t), with coefficients given by the moments 
at time t — 0: 



« oo ^ 

(121 



k=0 



Also, it is classically known (see Szego []) that both dctm„ and Pn{z) have integral 
representations; these facts lead to the following identities: 



rn{E,t) = f A{u)^l[pt{duk)=detmn{E,t) (13) 

k=i 

= detf(5,_,(t)) {pjk{E,0)) ) (14) 

\ ^ 0<i<n-l ^ 0<i<oo / 



0<j<oo 

0<j<oo 0<fc<n-l ' 



and 



Pn{E,t,z) = — / l[{z-Uk)A{ufY[pt{duk) (15) 



4 



using in the last equation 

eEr.»-fl_!^Ue^"(*-*)"". (17) 



z 



Equation (14) shows that the integral Tn{E,t) is the determinant of the projection of 
an infinite-dimensional plane onto the positive Fourier modes if+ = span{l, 2;, 2;^, ...}, 
as follows: 



n-l 



eEr*--' l^spanj ^ 0), A; = 0, n - l| © j ^ (18) 

Then, Sato [] tells us that this projection satisfies Hirota bilinear identities (Pliicker 
relations) and, in particular, the KP-hierarchy; but according to (13) this projection 
coincides with the integral Tn{E,t). Therefore Tn{E,t) satisfies for all n > the KP- 
equationsQ, 

i.j>0 

1.3 The vector (r„(_E, t))^>Q satisfies the Toda lattice 

Incidentally, equation (16) provides an elementary "classical" proof of Sato's repre- 
sentation of a wave-function in terms of a r-function. Indeed, in view of equation (23) 
below, we have that the L^— norm ^/hk of p, as in (10), has the form = rfc+i/r/j. 



so that pk = Pk\jhk = Pk\l T~k+i/Tk] this leads to Sato's representation of the so-called 
"wave vector" in terms of the "r-vector" r(t) := (r„(i?, t))„>o, 

^{E,t,z):={^k)k>o ■■= {pk{E,t,z))k>o 



e-2^T^^^'S{t)x{z), (20) 



A;>0 



where x(^) '■= (^*)«>oi this last equation holds, because Pkit, z) is a triangular linear 
combination of monomials z^ , represented by the (invertible) lower-triangular matrix 
S{t). Note zx = Ax, where A is the shift matrix A := (5i,j-i)i,j>o; therefore 

z"^ = L(t)\E', with L{t) = S{t)AS{t)^'^ symmetric, tridiagonal matrix. (21) 



'^Si{±d) = 5i (^±^, ±i^, ±i^, ...j , where Si{t) are the elementary Schur polynomials: 

00 
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The orthonormality of the polynomials pk{E,t,z) implies, by (10) and (20), the or- 
thonormality of the ^fc(-E, t, -2)'s for the weight po{dz) and thus taking partials in ti, 
one finds: 

— [ ^k{E,t,z)^e{E,t,z)po{dz)^0. 

oti Je 

Combining this relation with (20) and (21), one finds that ^' and L evolve according 
to^: 

^ = ^ins^ and ^ = l[{ns,L]. (Toda lattice) (22) 

For a full account, sec []. The seeds for these facts were already present in the 
pioneering work of Bessis, Itzykson and Zuber [] and later in the work of Witten []. 

1.4 The n-point correlation function 

The methods in this subsection are adapted from Mehta []. Return now to probability 
(8); upon expanding the products of the two determinants in the integral below, the 
denominator of (8) can be expressed as follows, using the orthogonahty of the monic 
orthogonal polynomials pk — Pk{F, t, z) on the full range F: 

f ^\z)f{pt{dz,) 

. n 

= / det(pi_i(2;j))i<ij<„det(pfe_i(2;^))i<fc,^<„ JJ Pt(d2;„) 

n=l 

n . 

= X] (-1)'''^''' n / PAk)-'i-i^k)P'K'{k)-l{Zk)Pt{dZk) 

7r,7r'e(Tra k=l ^ 

n—1 „ n— 1 

= nin / pl{z)pt{dz)^n\Jlhu. (23) 



Then, using the obvious fact (det A)"^ — det(A^'''), upon borrowing the exponen- 
tials from pt{dzi) and the hi from the denominator, one computes the probability (8) 
in terms of the polynomials p{z) := p{F, t, z) and a kernel Kn{y, z): 

P (spectrum M„ e £;") 

" ^im?. — L ^^M 5Z p3-i{zk)Pj-i{ze)] f[pt{dzi) 



^with regard to the Lie algebra splitting of g£(oo) into the algebras of skew-symmetric and 
lower triangular (including the diagonal) matrices Ai, (Borcl matrices): 

5^(00) = V,®Vb = A = A, + Ab. 
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If 

= — / det{Kn{zk, zi,))i<k,i<nT\po{dzi), (24) 
nl Je'^ ~ ~ 

where the kernel is defined in terms of the wave function (20) for E = F: 

n 

Kniy,z) ■.= J2^,.liF,t,y)^,_^{F,t,z). (25) 
The orthogonahty relations / '^kiF, t, z)'^i{F, t, z)po{dz) = 5ki lead to the reproduc- 

JF 

ing property for the kernel Kn{y, z): 

/ Kn{y,z)Kn{z,u)po{dz) = Kn{y,u), / Kn{z,z)po{dz) =n. (26) 
Jf Jf 

k 

Upon replacing E'^ by JJ^fi^j x F^''^ in (24), upon integrating out all the remaining 

1 

variables z^+i, ■■■■.Zn and using the reproducing property (26), one finds the n-point 
correlation function 

P(one eigenvalue in each [zj, Zi + dzi\, i = 1, ...,k) 

k 

= c„det {Kn{zi,Zj))^^.j^j^Y[po{dzi). (27) 

1 

Finally, by Poincare's formula for the probability P (^j E^ , the probability that no 
spectral point of M belongs to E is given by a Fredholm determinant 

oo „ k 

det(/ -\A) = l + Y.i->^)' / det(Aiz.,z,)) .^Mpo{dz.), 

k=i Jzi<...<Zk ^ ' ^ 

for the kernel A = K^, given by (25): 

P(no spectrum M e E) = det(/ - K^), K^{y, z) = K^^iy, z)Ie{z). (28) 

1.5 Toda vertex operators and boundary-time Virasoro constraints 

This subsection is taken from []; as a new ingredient, we introduced in [] the "vector 
vertex operator" for the Toda lattice, 

y.{t,z) ■=A-\{z^)e^i hz^^r^l^i — ai"^ (29) 

acting on the vector (r„)„>o; it can be viewed as a Darboux transformation, mapping 
a Toda solution into a new one. It satisfies the following simple "operator identity''^ 

+ '■])/^^X(t,z) = 0, £>-l, (30) 

4(J,XT)fc = (J,)fc(XT)fc, with (XT)fe = z2(fc-i)eEr*»^'e-'Sr ^^rfc_i 
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where the vector operator^ 

= ^ + ^(-A;)t„, + 4,ociiag(...,-2,-l,0,l,2,...), /c G Z 
J?^ = E:jf^jf^ keZ (31) 

i+j=k 

form the generators of the Heisenberg (oscillator) and Virasoro algebras respectively. 

The Hirota bilinear identities for r-functions, already mentioned in section 1.2, 
also generate so-called Fay identities, familiar to algebraic geometers; they are quadra- 
tic relations between translates of r-functions. These Fay identities enable us to ex- 
press the kernel Kn-, defined by (25), and the n-point correlation functions 
det(/r„(zi, Zj)), in terms of vertex operators acting on the r-vector r(t) := (t„(F, t))n>o, 
underlying \I/(F, t, z) as in (20). Namely, 

ir„(z,z) = X:^,-i(t,^)' = r-i(X(t,z)r(t))„, (32) 
i=i 

and, more generally, using higher degree Fay identities: 

det(if„(zi, 2;,))i<ij<fc = (\{ X(t, Zi)T{t)\ . (33) 

Vi=l / n 

Using the Neumann series in A for det(/ — XK^^), and setting E'^ := F\E, it follows 
that 

, TniE^,t) (e-/^^(*'^)^""'^''^V(t))„ , , 

P no spectrum e E) = \ , = 34 

Given the disjoint union 

r 

E= [j[A2i-i,A2i\cF CR, 

i=l 

the vectors r(t) := T{F,t) := (t„(F, t))„>o and T{E'^,t) := (r„(i?'^, t))„>o obey the 
following Virasoro-like constraints for m = —1, 0, 1, 



(a. jS+I - 6. jSiUi) r{F, t) = 0. (35) 



j>0 



- E + § (^^-^S- - &.JS1+i) I riE^^ t) = 0. (36) 



^the agreement is the following: — — = for fc < and tk = for fc > 0. Normal ordering 

otk 

means: pull the differentiation to the right, regardless of commutation relations. 
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A sketch of the proof of equations (35) and (36) goes as follows: besides the matrix 
L satisfying = z^^, we introduce a new matrix M, a matrix analogue of the 
Orlov-Schulman operator for KP, such that M\E' = {d/dz)'$; of course, [i^, M] = 1. 
The boundary condition (1) for the weight p{z)dz on F implies the vanishing of the 
integral 

Working out this integral and using the operators L and M, one proves the skew- 
symmetry of the matrix Q := Mf{L) + |(/' - g){L), and hence olQ [Q,L"'+^]^; 
remember / and g are the functions appearing in the weight (1); incidentally, the 
skew-symmetric matrix Q actually satisfies the string equation [L,Q] = f{L). We 
then use the ASV-correspondence [] between symmetry vector fields on the Toda wave 
functions and the Virasoro symmetries on the r-functions; this establishes equation 
(35). Equation (36) follows from (35), using (30) and taking into account the bound- 
ary of E"^. Observe the decoupling of the Virasoro constraints in (36) into a time-part 
and a boundary-part ! 

1.6 "Non-commutative" KP hierarchy for the probability distribution 

Let now the rational function V'{z), appearing in the weight (1), be of the following 
form: 

= I = ^i±^; (38) 

/ ao + aiz + a2Z'' 

then the equations (36) and their powers enable us to extract all partial derivatives 
Q'^i+'-'+^iji^E^ ,t) / dt\^ ...dt'l^ evaluated at t = in terms of partial differential opera- 
tors in the boundary points of the disjoint union E, 

A = E4'-'"/(^.)tti- A; = 1,2, 3,..., (39) 

where m = max(deg/ — 1, degg). Declare Ak to be of homogeneous weight k. 

When 02 = in V, the probability Vn{Ai, ...,A2r) '■= P(no spectrum M„ G E) 
satisfies the KP-hierarchy (19), but with the (commutative) partial derivatives in t 
replaced by the (non-commutative) partial differential operators Ak defined by (39)[]: 

Vn ■ AkAlVn - AkVn ■ AlVn - ^ Si{A)Vn ' S,{-A)Vn 

-|- (terms of lower weight i for 1 < i < A;) = 0, for A; > 3. (40) 
When a2 7^ in V, the ^^'s in (36) are replaced by more complicated expressions. 
6[^,S]t = ^{AB + BA) 

'''Si{zkA) :— Si{±Ai,±-A2,±-A3,...) are the elementary Schur polynomials with tk replaced 
1 2 3 

by -Ak- 
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Examples: When V{z) = z"^, the pfc(R, 0,^)'s are Hermite polynomials and, for a 
semi- infinite interval E = (A, oo), the first equation of the system (40) reduces to 
the Painleve IV equation for 9/9^4 log det(/ — K^). The case V — z — alogz leads 
to Laguerre polynomials at t = and to Painleve V for {A d/dA) logdct(J — K^). 
The Jacobi polynomials are also a special case of (38) and the first equation of the 
hierarchy (40) reduces to an unknown equation. These Painleve equations were first 
obtained by Tracy and Widom [] and then by us [], as a special instance of the 
non-commutative KP. 



2 Continuous kernels, KP and Virasoro 

Here, we shall be dealing with continuous kernels, like the Sine, Airy or Bessel kernels 
listed in (7). How do we obtain PDE's for the spectral distribution, using the theory 
of integrable systems ? The question is how to introduce time ? This question was 
resolved by Adler, Shiota and the author in []. In the finite situation, the kernel Kj^ 
in (25) was given by a Christoffel-Darboux sum, based on eigenvectors of L: 

n 

Kniy,z) ■.= J2'^,.iit,y)^j.i{t,z), (41) 

with := (^j)j>o an eigenvector of a second order problem and evolving in time 
according to the Toda lattice (see (21) and (22)): 

The continuous analogue goes as follows: consider now wave and adjoint wave func- 
tions '^{x,t,z) and '^*{x,t,z), with x e R, t e C~, z e C satisfying the KP- 
hierarchy. 



5* 

z^ = L*, — = (L")+^r, 



9^* 



Br 

— =[{L-U,L], (43) 



where L is a pseudo-differential operator 

d 

L = D + a_iD'^ + with D = —. 

ox 

We consider the p-reduced KP hierarchy, i.e., the reduction to L's such that — 
+ . . . is a differential operator for some p > 2. The precise continuous analogue 
of (41), which relates to a second order problem, would be to choose p — 2. 
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For the time being, we take p > 2 arbitrary. According to Sato's theory, we have 
the following representation of and ^P* in terms of a r-function (see []) 

nx,t,z) = e-+Sr^^-' ^(^-/f']) , and ^*{x,t,z) = e'^^-^T + l^-']) 

(44) 

for notation [z see (16). Observe^ 

^{x,t,z) := hu>^{x,t,ujz) and ^*{x,t,z) = ^ aa;\&*(x, t, cuz), (45) 

are the most general solution of the spectral problems = and L^p$* = 2;*'$* 
respectively. 

The continuous analogue of (41) is 

/X 
dx $*(a;, t, y)^{x, t, z), also K^{y, z) := K^{y, z)Ie{z), 

where $ and $* are given by (45), subjected to the condition X^^GCp ^77^ = 0- Intro- 
duce the p-reduced vertex operator 



z 



which maps the space of p-reduced KP r-functions into itself, with 

X{t,y,z) := ^ei:r{^-^')*»e^^(^"'-^"')^^; 
z-y 

see []. We now set c^^i = a^b^/ in the operator Y above; again using Fay identities, 
we obtain continuous analogues of (32), (33) and (34), 

K,{y,z):=-^Y{t,y,z)r{t) 



1 

det(^K^{yi,Zj))^^, .^^ = -]]_Y{t,yi,Zi)T 

det(/-Airf) = -e-^^'^^^(*'^'^V. 
r 

If a r-function r(t) satisfies a Virasoro constraint^] (analogue of (35)): 



r = CkpT for a fixed A; > — 1 



^(p := {lj such that ujP ~ 1} 

""4'^ = J; + (-^)*-^ JP = J2^+J=i ■■ J't^jf^ see footnote 4. 
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then, given a subset E — [J [A2i_i, ^21] C R+, the Fredholm determinant of 



^"^(A, A') := -%|^/^(A'), A = z^, X' = z^ 
P z 2 z 2 

satisfies the following constraint 

/ 2r 8 1 \ 

-p E + 2 ^-^S^ " ''^^ ^ ^'^^^ " ^^^^^ ^ °- 



1=1 



Notice again the decoupling of the boundary- and the time-parts, as in (36). We 
further proceed in the same way as in the discrete case: Under suitable assumptions 
on the initial conditions, analogous to (38), we can express all t-partials at t = in 
terms of partial differential operators in the endpoints Ai oi E, provided p = 2; we 
then substitute those expressions in the KP hierarchy, thus leading to a hierarchy of 
PDE's in terms of non-commutative operators in the Ai. A precise statement due to 
Adler-Shiota-van Moerbeke [] proceeds as follows: 

Theorem. Consider a Grst-order differential operator A in z of tlie form 

A = A, = \z--^' f I- + Viz)] + E c,z-'\ (46) 



witii 



9z , ,>i 



V[^z) = -z + -z^^O, m = deg V = or 2, (47) 
2 6 



and its "Fourier" transform 

A^A,^(^^{x + V'{D)) +1 + 5^ j witli D^-^. (48) 

Let ^{x, z), X E 'R, z E C be a solution of the linear partial differential equation 

A,^{x,z) ^ A^<lf{x,z), (49) 

with holomorphic (in z~^) initial condition at x — 0, subjected to the following 
differential equation for some a,b,c E C, 

{aA^ + bA + c)*(0, z) = ^^^(0, z), with *(0, z)^l + 0{z-^). (50) 

Then 

• ^(x, z) is a solution of a second order problem for some potential q{x) 

(D^ + q{x))<i/{x,z) ^ z'^<if{x,z). (51) 
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• Given a subset E = [J [^21-1, A2i] and the kernel 

i=l 

K^{y,z) := Ie[z) j ^TTZ^TTi (^^j 

with 

w=±i 

tiie Fredholni determinant V{Ai, A2r) '■= det(/ — XK^) satisfies a hierarchy of 
bihnear partial differential equation^ in the Ai for odd n > 3: 

V-AnA{P-AnV-A{P- J2 S,{A)r ■ Sj{-A)r 

i+j=n+l 

+ (terms of lower weight i for I < i < n) = 0, (53) 
where the An are differential operators of homogeneous "weight" n, defined by 

^n = E^. " TTT^ n = l,3,5,.... 

i=l '^^i 

examples. Airy kernel: Given an entire function U{u) growing sufficiently fast at 
cxD, consider its "Fourier" transform, 

/CXD 
-00 

define the associated function p{z) and the differential operator A: 

p{z) = -i=e^(^)-^'(^)v^ A := p(^)^^pW~'- (55) 



Then the function 

^{x, z) := p{z)F{x + U'{z)) (56) 

obeys, for a trivial reason, 

d 

A^(x,z) = — ^(x,z), (57) 
ox 

and, upon integration of ^e"^^"-*"^"^, one shows '^{x,z) also satisfies: 

(u'i-^) - x\ ^(x, z) = U'{z)^{x, z). (58) 



io5,(±^) :=5,(±A,0,±i^3,0,...) 



13 



Upon letting x — in (58) and using (57), we find 

U'{A)^{0, z) = f/'(z)*(0, z). (59) 

Choosing U{z) = z^/3, the integral (54) becomes the Airy function; the operator A 
in (55) coincides with (46), for V{z) = ^z^, a = 0, /? = 4, m = 2: 

A:= — (^ + 2zA - -z-^ and A = D. 
2z \dz I 4 

Moreover (59) becomes equation (50) and, by stationary phase, satisfies the asymp- 
totics requested in (50): 

^2^(0, z) = ^2^(0, z), with ^(0, z) = l + 0{z-^) for 2/00. 

From (58), it follows that q{x) = —x in (51); the kernel (52) reduces to the Airy 
kernel mentioned in (7) and therefore its Fredholm determinant satisfies hierarchy 

n-l 

(53) with An = Sill Ai ^ d/dAi. When is a semi-infinite interval, the first PDE 
in the hierarchy (53) reduces to Painleve II for {8/ OA) logdet(J — K^), as originally 
discovered in []. Note that in this example, the r-function is the Kontsevich integral. 
For details on our methods, see []. 

Bessel kernel: Pick ■m^V{z)^ A^ and A^ in (46), (47) and (48) as follows, 

V{z) = -z, with m = 0, A := A, := ^2;(^ - 1), i := := ^(x - 1)^, 
and consider the partial differential equation (49) with initial condition (50): 

AA^ - 2A - + ^(0, z) = ^2^(0, z), with ^(0, z) = 1 + 0(2"^). 

(Bessel equation) 

The solution to this problem is given byQ 
and 

^{x,z) = e'^'Biil- x)z). (60) 
The latter is an eigenfunction for the second order problem 

(d'-^^^]^{x,z) = z'^{x,z). 
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Choosing in (52), b+ = e~^/2y^, 6_ = ib+, and integrating from 1 to a; = i + 1 
leads to the Bessel kernel mentioned in (7). If £■ is a semi- infinite interval, the first 
equation of (53) reduces to the Painleve V equation for (—Ad/dA) logdet(/ — K^), 
which was first done in []. Specializing to z/ = ±| leads to the sine kernels, mentioned 
in (7). Finally, the r-function associated with the wave function is a double Laplace- 
like transform; see []. The details for the Bessel kernel can be found in []. 
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